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Certain regularities of the motion of glass melt flows in a highly efficient furnace for container glass are stu¬ 
died employing the method of physical modeling of the glass melt hydrodynamics. 


Modern practice calls for studying the conditions that en¬ 
sure production of high-quality glass in highly efficient 
glass-melting furnaces. 

The method of physical modeling [1] allows for a clear 
picture of mass exchange in a glass-melting furnace and 
makes it possible to issue recommendations for furnace ope¬ 
ration. The present paper describes the glass melt hydrody¬ 
namics in the melting tank of the furnace located in the 
Krasnoe Ekho glass factory. 

Technical parameters of the furnace 

Type of furnace.Regenerative, with transverse 

flame direction and 
with working channels 

Glass.Clear glass BT-1 

Furnace output, ton/day.200 

Surface area, m 2 : 

of the melting tank.104.93 

of the working channels.18.78 

Type of fuel.Natural gas 

(gcak = 33 MJ / m 3 Qr 7970 kcal / m 3 ) 

Product range.Bottles for various purposes 

Type of glass-shaping machines . . AL-118-2 (3 machines) 

Surface area of entries, m 2 .5.86 

Fuel consumption for furnace heating, nm 3 /h . . . 1801.55 
Temperature conditions, °C: 

max. melting temperature.1530 

max. working temperature.1220 

air heating temperature in regenerator.1100 

The rectangular melting tank has a trapezoidal spillover 
threshold and a sunk clarification and homogenization zone 
(a deep refiner). 

The furnace is equipped with two working channels with 
independent heating. The melting-tank bottom, including the 
deep refiner, is heat-insulated by two chamotte layers whose 
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total thickness is 600 mm. The melting-tank walls have con¬ 
tinuous chamotte heat insulation. The furnace is heated by 
four pairs of gas burners. The fuel is supplied to the burners 
via gas nozzles with an exit diameter of 120 mm, which are 
arranged facing each other at an angle of 90°. 

The physical modeling method is based on similarity the¬ 
ory, according to which dimensional physical values can be 
combined in dimensionless complexes in such a way that the 
number of complexes would be smaller than the number of 
values composing these complexes [1,2]. Furthermore, simi¬ 
larity theory determines the conditions under which the mo¬ 
deling results can be extrapolated to other processes that are 
similar to the ones studied. 

The calculation of the principal parameters for the mo¬ 
deling was carried out according to the method developed by 
A. A. Sokolov and K. A. Pchelyakov [3]. 

In accordance with the problem, the following tempera¬ 
ture range was chosen for the glass melt: £ max = 1500°C, 

^m m = 1370 ° C , Af= 130°C. 

Figure 1 shows the temperature dependence of the kine¬ 
matic viscosity of the glass melt and the model liquid for the 
glass melted in glass-melting furnace No. 2 and the tempera¬ 
ture relation plot. 

We find from the plot that the maximum temperature 
correlates with the minimum kinematic-viscosity value 
0.32 x 10“ 3 m 2 /sec. 

The geometric scale is related to the viscosity scale. 
Based on the following main conditions of the modeling of 
glass melt convection: 

Gr - gl 3 /v 2 - idem; 

(3/A/ = idem; 

Re = v//v = idem, 

we obtain 

c — sir 2 • c =r 3 / 2 

W V V 5 L V ’ 
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Fig. 1. Temperature-viscosity parameters of glass melt (a) and 
model liquid (b ); temperature relation plot ( c ). 

where Gr is the Grashof number; g is the free-fall accelera¬ 
tion, m/sec 2 ; / is the geometric size, m; v is the kinematic 
viscosity, m 2 /sec; P is the volume expansion coefficient, 
K _1 ; Re is the Reynolds number; v is the velocity, m/sec; 
C is the scale. 

The scale C } = 25 was chosen for modeling, and then 
C v = 125 (Cj is the linear scale, C v is the viscosity scale). 
Taking these scales into account, we will find the maximum 
and minimum values of the kinematic viscosity of the model 
liquid: 

v m = 0.18 x 10 -3 m 2 /sec; v m = 0.026 x 10 -3 m 2 /sec. 

The indicated kinematic-viscosity values correlate with 
the following temperatures of the model liquid: t ™ m = 48°C, 

t m = 85°C, A/ m = 37°C. The temperature scale: C, = 

At/At m = 130/37 = 3.5. 

The scale of velocity is found based on the condition 
Re = idem: C v = C v /C [ = 5, and the output scale is found 
from the condition C G = C p C v C l = 5937 (C is the density 
scale). The consumption of the model liquid amounts to 

G m = G/G g = 23 g/min, 

where G is the furnace output, ton/min. 



Fig. 2. Experimental device for modeling of glass melt hydrody¬ 
namics. 

The time scale C T = = 5. Glycerin was used as the 

model liquid, with LiCl used as an additive. 

The model design has to ensure not only the similarity of 
the thermal and hydrodynamic processes corresponding to 
the processes in the functioning furnace but also convenience 
of observation and the ability to measure the main flow pa¬ 
rameters: glass melt routes and speed and temperature fields. 

The scheme of the model is shown in Fig. 2. The case 1 
of the model is made of organic glass 4 mm thick, and the 
roof 2 is of sheet aluminum. The model is mounted on a spe¬ 
cial stand. Water-cooled slabs 3, 4 positioned beneath the 
bottom of the model are connected to thermostats and corre¬ 
spond to two furnace zones: melting zone and homogeniza¬ 
tion/clarification zone. The temperature in the thermostats 
was maintained with an accuracy to 0.5°C. By varying the 
temperature in the thermostats and the consumption of ther¬ 
mostat liquid, it is possible to regulate the temperature in 
each zone. 

Heating of the model liquid from above was performed 
by Nichrome heaters 5, which made it possible to control the 
heat release level within a wide range. Each Nichrome heater 
section was connected to an autotransformer. The heat flow 
was controlled by varying the voltage applied to the heater, 
and that is how the temperature distribution on the model liq¬ 
uid surface was controlled. 

Glass production was simulated by pumping the model 
liquid from the working tank 6 using the gear pump 7. The 
liquid was cooled in the transport pipeline and the chiller 8 to 
the temperature of the glass melt in the zone below the batch. 
Thermocouples were connected to the digital voltage meter 
9. The total absolute error of a thermocouple and the digital 
meter did not exceed 0.5°C, and the relative error amounted 
to 1%. 

The hydrodynamic processes in the furnace model were 
studied by registering the migration routes of colored indica¬ 
tors introduced into the liquid in the form of a vertical fila¬ 
ment employing a special syringe. 
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The effect of the following factors on the glass melt hy¬ 
drodynamics was studied in the first stage: position of the 
maximum-temperature zone (quellpunkt) on the glass melt 
surface, specific output of glass melt, and height of the 
spillover threshold. 

Specific efficiency (glass melt output) of the furnace. 

Model studies were conducted within the productivity range 
of 1500 - 3000 kg/m 2 per day. 

When the specific output is equal to 1500 kg/m 2 per day, 
in addition to the main flows, a circulating flow known as the 
Schuld wave exists in the zone following the spillover 
threshold, and another flow exists in the thin surface layer. 
These latter flows are regarded as undesirable and capable of 
degrading the thermal homogeneity of the melt. When the 
output is higher than 1500 kg/m 2 per day, such flows are not 
registered. An increase in the glass melt output from 2000 to 
3000 kg/m 2 per day did not significantly alter the melt hy¬ 
drodynamics either in the melting zone of the furnace or in 
the clarification zone. The speed of the glass flow in the bot¬ 
tom layer in the melting zone and the speed of the flows ad¬ 
jacent to the walls in the clarification zone increased (Fig. 3). 

Probably, with the specific output exceeding 2000 kg/m 2 
per day, it is necessary to provide for a higher concentration 
of thermal energy in the spillover threshold zone. This can be 
accomplished, for instance, by installing additional electric 
heating in front of the spillover threshold and directly behind 
this threshold (to decelerate the glass melt flow that migrates 
along the wall from the threshold to the neck). 

Quellpunkt position. The main parameters of the fur¬ 
nace performance: glass melt output 2000 kg/m 2 per day, the 
tank bottom is heat-insulated, melting temperature 1530°C. 

The migration of the glass melt in the tank is indicated in 
Fig. 4 for different positions of the quellpunkt. An analysis of 
these data first of all reveals a significant effect of the posi¬ 
tion of the maximum-temperature zone on the glass melt mi¬ 
gration. Obviously, special attention should be paid to this 
factor when searching for the optimum operating conditions 
for the furnace. 

The best option among the three operating modes is the 
one in which the quellpunkt is located above the threshold or 
in the direct vicinity of it. In this case, the pouring cycle of 
the glass melt flows is longer, and the melt flows at the en¬ 
trance to the clarification zones are more stable. 

Threshold height. The height of the spillover threshold 
in the functioning furnace is 270 mm. The modeling demon¬ 
strated that with this threshold height, the reverse current 
flowing underneath the batch has a small thickness and insuf¬ 
ficient speed. When the threshold height was increased to 
600 - 700 mm, the reverse-flow thickness amounted to about 
1/2 of the melting-tank depth. At the same time, the speed of 
the flow increased 1.4- 1.6 times, which ultimately should 
have a positive effect on the glass melting process. 

The second phase of modeling involved determination of 
the optimum operating mode for a glass-melting furnace 


Quellpunkt 



Fig. 3. Glass melt currents in model of glass-melting furnace of ca¬ 
pacity 2000 kg/m 2 per day (a) and 3000 kg/m 2 per day ( b ). Num¬ 
bers at the arrows: speed of the flows (m/h). 


Quellpunkt 



Fig. 4. Glass melt flows in model of furnace with various positions 
of the quellpunkt. 

whose melting tank was divided by the threshold into two 
parts: the melting and clarification zones. 

The optimum operating conditions were chosen based on 
the following assumptions [4]: 

- intense circulatory movement should take place in the 
melting zone; 

- all or a major part of the glass melt should traverse the 
maximum-temperature zone; 

- a stable and clearly expressed quellpunkt should exist 
in the melting zone; 

- the speed of the upper branch of the pouring cycle 
should be as high as possible. 

The experiments showed that the qualitative aspect of the 
glass melt currents in the furnace significantly depends only 
on the degree of heat insulation of the furnace bottom. When 
the bottom is not insulated and the temperature difference 
across the depth is 120 - 140°C, two circulating flows coex¬ 
ist in the furnace: the pouring flow and the working flow, 
whose respective thickness amounts to 2/3 and 1/3 of the 
melting zone depth. Besides, another working current exists 
near the bottom and is directed toward the charging zone. A 
low-mobility glass melt layer exists between the pouring and 
working flows. 

Depending on the location of the temperature maximum, 
the extension of the flows can vary, but on the whole the hy¬ 
drodynamic picture is preserved. In this way, three flows ex¬ 
ist in the melting zone: the upper flow (it is directed under¬ 
neath the nonmelted batch and maintains the boundary be¬ 
tween the batch and the froth), the middle flow (the lower 
branch of the pouring cycle), and the bottom flow. The mid¬ 
dle flow rises in the quellpunkt area and splits into two 
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branches: one of them completes the pouring cycle, and the 
second branch moves inside the surface layer toward the tank 
neck. The bottom flow descends along the vertical wall of 
the clarification zone and arrives at the working zone. 

In the case where the furnace bottom is heat-insulated 
and the temperature difference across the glass melt depth is 
about 70 - 80°C, the hydrodynamic situation in the furnace 
is different. Only one joint pouring cycle exists in the melt¬ 
ing zone, the upper flow of which is directed toward the 
charging zone, and the lower one toward the working zone. 
The lower flow rises in the quellpunkt zone and splits into 
two branches, one of which completes the pouring cycle, and 
the other moves toward the clarification zone. The speed of 
the bottom glass layers increases 2.5-3 times in the case of 
heat insulation of the bottom, which will make it possible to 
intensify the melting process. 

The study performed resulted in the following conclu¬ 
sions: 

- the spillover threshold contributes to increasing the 
clarification potential of the furnace and decreases the 
amount of low-mobility glass melt; 

- the optimum thermal conditions in this case are accom¬ 
plished if the quellpunkt is located directly at the threshold or 
above the threshold; 

- the trapezoidal shape of the threshold facilitates the 
thermal homogeneity of the glass melt and extends the 
threshold service life; 


- the described results of the modeling substantiated the 
need to insulate the melting zone of the furnace tank; this 
produces a two-level convection, which makes it possible to 
identify the quellpunkt with accuracy; 

- with the output exceeding 2000 kg/m 2 per day, the 

number of low-mobility melt sites decreases significantly; 

- the hottest glass melt currents circulate near the neck 

wall; 

- it is expedient to increase the threshold height to 
0.7 - 0.9 m for rational organization of the flow and intensi¬ 
fication of the reverse flow oriented underneath the batch. 
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